AO-AOM  213 


UNCLASSIFIED 


AIR  FORCE  INST  OF  TECH  WRI6HT-PATTERS0N  AFB  OHIO  SCH— ETC  F/6  20/5 
ANALYTIC  electron  ENERGY  DISTRIBUTION  FUNCTIONS  FOR  ELECTRIC  OI~ETC(U) 
DEC  78  T E BIST 


AFIT/GEP/PH/TBD-A 


NL 


A FI T/GEP/PH/78D- 4 


ANALYTIC  ELECTRON  ENERGY  DISTRIBUTION 
FUNCTIONS  FOR 
ELECTRIC  DISCHARGE  LASERS 


THESIS 


AFIT/GEP/FH/78D-4  Thomas  E.  Gist 

Capt  USAF 


Mpp  c wf  A pabl>c  Y f.  t-  e wif  J u.\  1 1 


.k.  Jk,  %J 


AFI T/GEP/PH/7  8 D- 4 


(J^ 


ANALYTIC  ELECTRON  ENERGY  DISTRIBUTION 

■ . r * 

FUNCTIONS  FOR 


ELECTRIC  DISCHARGE  LASERS, 


TiffiSIS - 

(^f  ^ r, 

Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science 


Thoraa s E y/ci stj  B.S. 


Capt 


USAF 


Graduate  Engineering  Physics 


»78 


7oCf>- 


Approved  for  public  release:  distribution  unlimited. 


Acknowledgemnts 


I would  like  to  acknowledge  the  guidance  and  assistance 
of  Philip  Nielsen,  who  provided  the  impetus  for  this  thesis, 
and  without  whose  help,  completion  would  have  been  very  dif- 
ficult. I would  also  like  to  acknowledge  the  advice  and  as- 
sistance of  Allen  Hunter,  who  on  several  occasions  provided 
the  insight  necessary  to  get  past  difficulties.  I would  like 
to  acknowledge  the  aid  and  support  of  William  Ercoline,  Harold 
Hastings,  and  Edward  Seward.  On  many  occasions  they  were  able 
to  provide  information  they  had  discovered  that  saved  me  much 
time  and  effort.  Finally,  I must  acknowledge  the  support, 
aid  and  understanding  of  my  wife,  Karen.  Her  long  hours  of 
typing  and  plotting  were  instrumental  in  the  completion  of 
this  thesis?  her  understanding  support  was  vital. 


Thomas  E.  Gist 


Contents 


Acknowledgemnts  

List  of  Figures  

List  of  Tables  

Abstract . 

I.  Introduction  

Background  

Problem  Statement  

Definitions  

Assumptions  

Scope 

Summary  of  Current  Knowledge  

Approach  

II.  Solution  of  the  Time  Independent  Boltzmann 
Equation  for  Inelastic  Collisions  . . . . 

Rationale  for  Approach  

Derivation  of  the  Equation  for  T(€)  . . . 
Boltzmann  Equation  for  F(€)  

Form  for  € and  R(€)  

Derivation  of  Equation  for  T(€)  . . . 
Solution  for  f(6),  Two-Level  Case  . . . . 

e>x 

€<x 

Complete  Solution  

Solution  for  f(€),  Three-Level  Case  . . . 

e>x 


Complete  Three-Level  Solution 
Extension  of  the  Solution  .... 
More  Levels  (X^^  < 

Xj^  > (X2  - 

III.  Analysis  and  Conclusion  

General  Considerations  

Criteria  for  Evaluation  

Theoretical  Cross-Sections  . . . 

Two-Level  Model  

Three-Level  Model  

Realistic  Cross-Sections  . . . . 

Two-Level  Model  

Three-Level  Model  

Modeling  a Multi-Level  System 
Superelastic  Collisions  . . . 


iii 


General  Considerations  . . 

Physical  

Mathematical  

Two-Level  Model  

dN/dq  =0  

Continuity  of  df/d€  at  q 
Validity  of  Assumptions 
Three-Level  Model  


Appendix  C:  NGB  

General  Description  . . 
Particular  Application 
Program  Modification 
Input  Parameters 

Appendix  Ds  TI-59  Programs  . 

General  

Two-Level  Program  . . . 
Program  Listing  . . 
Program  Use  . . . . 
Three-Level  Program  . . 
Program  Listing  . . 
Program  Use  . . . . 


iv 


m ‘ 

- 

List  of  Figures 

■j 

i 

Figure 

Page  i 

1 

Published  Distribution  Functions  in  Ng  . . . 

10 

2 

Analytic  R and  

12 

3 

Comparison  of  threshold  energy  relationships 

22 

4 

Ideal  Two  Level  System,  X/kT  =10  

27 

5 

Ideal  Two  Level  System,  X/kT  =20  

X 

28 

6 

Ideal  Two  Level  System,  X/kT  =30  

X 

29 

7 

Ideal  Two  Level  System,  X/kT  = 100  

X 

30 

8 

Demonstration  of  Finite  Differencing 

Discrepancies  

32 

9 

Ideal  Three  Level  System  

33 

10 

Two  Level  System,  El  Level  of  N2  

36 

11 

Two-Level  System,  E2  Level  of  

37  ; 

12 

Two  Level  System,  E2  Level  in  Argon  

00 

13 

Three  Level  System,  El  and  E4  Levels  of  N2 

40  1 

14 

Three  Level  System,  E2  and  E4  Levels  of  N2 

41  1 

15 

Three  Level  System,  El  and  E5  Levels  of  N2 

42  1 

16 

Five  Level  System,  V6,  V7.  El,  and  E2  Levels 

in  CO 

44 

17 

Effect  of  Superelastic  Collisions  

46 

B-1 

Comparison  off(€)  €^0  and  € X . . . . 

62 

B-2 

Comparison  of  f(G)  € » X^^  and  € « X2  • • • 

65 

1 

L 

V 

t 

.i 

. m piHiiii' m hi  iimii  ii  ii 

List  of  Tables 

Table 

Page 

B-1 

Convergence  of  3f/a€  

D-1 

Program  for  Two-Level  Model  .... 

68 

D-2 

Program  for  Three-Level  Model  . . . 

74 

D-3 

Initial  Parameters  

88 

D-4 

Post-Initialization  Parameters  . . . 

88 

AF I t/GEP/PH/7  8D - 4 


Abstract 

The  Boltzmann  equation  for  electrons  is  solved  for  back- 
ground gases  of  one  or  two  excited  levels  under  conditions 
typical  of  electric  discharge  lasers.  Reaction  rates  are 
assiuned  constant  above  a threshold  energy,  and  the  electron 
heating  rate  is  assumed  independent  of  electron  energy. 

Elastic  and  superelastic  collisions  were  not  considered.  Nu- 
merical solutions  from  a program  developed  in  conjunction  with 
Seward,  Hastings,  and  Ercoline  are  in  excellent  agreement  with 
the  analytic  solutions.  The  analytic  solutions  are  compared 
to  numerical  solutions  that  use  realistic  cross  sections:  the 
analytic  solutions  agree  well  for  some  cross  sections.  It  is 
possible  to  determine  a priori  which  realistic  cases  the  an- 
alytic model  will  fit.  For  these  cases  it  is  possible  to 
estimate  pumping  rates  for  upper  levels  with  reasonable  suc- 
cess. In  some  cases  it  is  possible  to  model  systems  with 
more  than  two  excited  levels.  Limitations  to  the  analytic 
model  and  validity  of  the  assumptions  are  discussed.  Pro- 
grams for  evaluating  the  analytic  solutions  on  a Texas  Instru- 
ments programmable  59  calculator  are  included. 
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ANALYTIC  ELECTRON  ENERGY  DISTRIBUTION  ^ 

FUNCTIONS  FOR 
ELECTRIC  DISCHARGE  LASERS 

I Introduction 

1 

{ 

Background  ‘ | 

There  is  a considerable  interest  in  the  United  States 
Air  Force  in  developing  high  energy  lasers.  One  of  the  prob- 
lems that  must  be  solved  in  this  development  process  is  pre- 
dicting the  electron  energy  distribution  function  produced  by 
some  heating  machanism  in  the  plasma  in  a laser  gain  tube. 

To  date  this  has  been  attempted  both  numerically  and  analyti-  ‘ | 

cally,  without  complete  success. 

The  numerical  methods  will  produce  a value  for  the  num- 
ber of  electrons  with  energies  between  € and  € + d€  at  a i' 

given  time,  for  a given  set  of  initial  conditions.  But  it  is 
then  impossible  to  determine  from  this  solution  the  same  in- 
formation  for  a different  set  of  initial  conditions.  In 

addition,  without  having  a second  independent  solution  to  the  | 

problem  for  at  least  some  cases  that  can  be  compared  to  the  I 

numerical  solution,  it  is  difficult  to  be  certain  that  the  | 

numerical  results  are  valid.  i| 

The  existing  analytic  solutions  also  lack  completeness  j 

for  several  reasons.  It  is  not  always  possible  to  obtain  an-  ? 

t 

alytic  solutions.  Those  solutions  that  do  exist  are  based  ' 

on  versions  of  the  Boltzmann  equation  that  have  been  simpli- 


fied  by  various  assumptions:  although  these  assumptions  make 
solving  the  equation  simpler,  they  also  reduce  the  accuracy 
of  the  solutions.  In  addition,  these  assumptions  make  the 
solutions  valid  only  for  particular  classes  of  initial  condi- 
tions . 

The  analytic  models  are  very  important,  however.  Al- 
though they  are  not  valid  for  all  physical  situations  due  to 
the  various  assumptions  involved,  there  will  be  physical  sit- 
uations of  interest  that  the  analytic  solutions  will  model 
accurately.  In  these  cases,  the  analytic  solutions  can  be 
used  to  verify  the  validity  of  the  numerical  solutions.  In 
addition,  the  analytic  models  can  indicate  how  the  numerical 
solutions  change  as  the  physical  situation  changes.  The  an- 
alytic solutions  can  also  give  very  important  insight  into 
the  physical  phenomena  involved. 


Problem  Statement 

Analytic  solutions  of  the  Boltzmann  equation  are  devel- 
oped which  describe  the  electron  distribution  function  in  a 
plasma.  The  solutions  developed  eliminate  some  of  the  prob- 
lems associated  with  previous  analytic  solutions  and  are  com- 


pared to  a series  of  numerical  solutions  of  both  idealized 


and  realistic  systems. 

Definitions 

Several  definitions  are  needed  for  this  discussion.  The 


first  of  these  is  the  electron  distribution  function  f(€). 
If  F(€)d€/N  is  the  fraction  of  electrons  with  energies  be- 


II 


2 


r 


This  form  of  the 


tween  € and  € + df , then  f(€)  = - 
distribution  function  is  chosen  to  simplify  some  of  the  out- 


put for  the  numerical  solutions.  The  normalization  factor  N 
simplifies  the  derivation  and  is  defined  so  that 


/®  i ^ CO 

Jo 


d€  = 1 


(1-1) 


The  appropriate  Boltzmann  equation  for  F(€)  is 


3F(g , t ) 

at  a€ 


2 ^/F(g,t)  _ g aF(g,t) 

3 I 2 ^ as 


E 

3 


nj3^F(€’ .t)R|^(e’ ) - ngF(€,t)Rp(€) 


(1-2) 


where 


€ = rate  of  energy  increase  per  electron  per  unit 
time(electron  heating  rate) 
n^g  = number  density  of  species 
n^*  = number  density  of  species  3*  which  is  "inverse" 
to  3 

Rp(6)  = rate  constant  for  the  reaction  e(€)  +3  

e(e'  ) + 3* 

Rg(€')  = rate  constant  for  the  reaction  e(t')  + 3*  ^ 

e(€)  + 3 

3 and  3*  3-^®  "the  initial  and  final  states  of  some  transi- 
tion. 

e(€)  is  an  electron  of  energy  € (Ref  13 j 2348-2349) 

Another  term  that  needs  defining  is  "steady-state."  For 
the  purposes  of  this  discussion,  a steady  state  solution  for 
the  distribution  function  is  one  for  which  the  distribution 


3 


function  is  a function  of  energy  only,  not  time,  over  time 
scales  of  interest.  Note  that  this  does  not  require  a con- 
stant electron  density,  but  only  that  the  proportion  of  elec- 
trons in  any  given  energy  range  is  constant. 

Assumptions 

Several  assumptions  are  used  in  this  thesis.  First,  it 
is  assumed  that  a steady  state  condition  exists  after  some 
initial  transient  phase;  this  phase  typically  lasts  for  a 
time  on  the  order  of  picoseconds.  Second,  ix  is  assumed  that 
the  electron  density  is  low  enough  (<  10  cm  that  electron- 
electron  interactions  can  be  ignored.  These  two  assumptions 
are  direct  consequences  of  the  limitations  in  scope  to  be 
discussed  later.  Third,  it  is  assumed  that  f(€)  is  of  the 
form 


f(e) 


e > Xr 


(1-3) 


where  c is  an  arbitrary  constant  and  k is  Boltzmann’s  con- 
stant. Note  that  this  differs  from  the  normal  Maxwellian 
distribution  in  that  instead  of  T being  a parameter  charac- 
teristic of  the  distribution  as  a whole,  it  is  a function  cf 
the  electron  energy  €.  It  is  further  assumed  that 


€ » kT(€)  , € > X 


e 


II  » kT^  , € > X, 


d€ 


3 


(1-4) 

(1-5) 
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The  last  three  assumptions  (Eqs  (1-3).  (1-^),  and  (1-5)), 
are  based  on  observations  of  published  distribution  functions, 
as  discussed  in  Chapter  II.  Finally,  it  is  assumed  that 


where  n^  is  the  neutral  particle  density 
n^  is  the  excited  particle  density 

This  assumption  is  required  to  enable  super- elastic  colli- 
sions to  be  ignored.  For  gases  at  atmospheric  pressure, 
n^  can  approach  10^*^  - 10^^  cm~^  before  this  assumption  be- 
comes invalid.  Although  these  values  can  be  exceeded  in 
I electric  discharge  lasers  in  some  circumstances,  the  assump- 

tion would  still  be  valid  until  n exceeded  those  levels. 

The  validity  of  Eqs  (1-4)  through  (1-6)  is  addressed  in 
more  detail  in  Appendix  A. 

Scope 

This  thesis  examines  the  distribution  of  an  electron 
gas  in  the  presence  of  a background  gas.  It  deals  with  the 
steady-state  case  only.  It  includes  background  gases 
; with  both  single  and  multiple  excited  levels.  Superelastic 

collisions  will  not  be  considered  analytically  except  to 
determine  the  limits  of  validity  of  the  solutions,  and  elec- 
tron-electron collisions  are  excluded. 

Summary  of  Current  Knowledge 

At  present  there  are  few  published  works  that  deal 

i 


1 


exactly  with  this  problem,  although  there  are  several  that 
deal  with  closely  related  problems.  For  example,  Rockwood 
did  considerable  work  in  developing  a numerical  method  of 
solving  the  Boltzmann  equation,  but  did  not  find  any  analytic 
distribution  functions  (Ref  13).  Nielsen,  Canavan,  and  Rock- 
wood  used  a quantum-kinetic  approach  to  investigate  gas 
breakdown  numerically  (Refs  10;  11;  and  14).  Englehardt 
and  Phelps  used  numerical  techniques  to  solve  the  Boltzmann 
equation  in  order  to  find  excitation  cross-sections  (Ref  4). 
Elliott  and  Greene  used  numerical  techniques  to  study  the 
time  evolution  of  the  Boltzmann  equation  (Ref  3)-  Nighan 
used  numerical  techniques  to  find  distribution  functions  in 
various  gases  (Ref  12).  Frost  and  Phelps  were  able  to  solve 
the  Boltzmann  equation  analytically  for  rotational  excitations, 
but  only  by  assuming  that  the  rotational  levels  were  so  close 
together  as  to  form  a continuum;  this  assumption  is  invalid 
for  vibrational  or  excitational  levels.  For  these  levels 
Frost  and  Phelps  used  numerical  techniques  (Ref  5)-  Long, 
Bailey,  and  Garscadden  explicitly  consider  a spatially  non- 
homogenous  problem  by  using  a first-order  Laguerre  expansion 
of  the  Boltzmann  equation.  They  then  use  the  zero-order 
function,  f ^ , (which  corresponds  to  f in  this  paper)  to  find 
fj^.  However,  since  their  ultimate  goal  was  drift  velocities, 
which  they  could  find  without  knowing  f^  explicitly,  they 
did  not  find  an  analytic  distribution  function  corresponding 
to  f (Ref  7). 

Nielsen  has  one  unpublished  paper  which  addresses  this 
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specific  problem.  In  it  he  uses  the  same  assumptions  as  this 
thesis  to  find  an  analytic  form  for  f(€),  suggests  how  the 
solution  might  be  extended,  and  discusses  the  impact  of  his 
solution  on  numerical  programs.  However,  this  thesis  shows 
there  is  a minor  error  in  Nielsen's  paper;  furthermore, 
Nielsen's  paper  does  not  actually  extend  the  solution  to 
three  or  more  levels,  nor  does  it  include  any  extensive  com- 
parisons with  numerical  solutions  (Ref  9)* 


Approach 

This  thesis  is  an  extension  of  Nielsen's  unpublished 
paper  (Ref  9)*  Using  the  same  general  method,  a corrected 
two  level  model  is  developed.  This  model  is  then  extended 
to  three  levels,  where  the  spacing  between  the  two  excited 
levels  is  greater  than  the  spacing  between  the  ground  level 
and  the  first  excited  level.  An  example  of  a real  system 
with  these  characteristics  would  be  one  with  a low-lying 
group  of  closely-spaced  vibrational  levels,  which  could  be 
approximated  by  a single  lower  excited  level,  and  then  a 
higher  electronic  level.  Then  it  is  shown  how  the  solution 
can  be  extended  to  more  than  three  levels,  with  arbitrary 
spacing. 

While  the  analytic  research  was  being  done,  a numerical 
program  was  also  being  developed  in  conjunction  with  E.  D. 
Seward,  W.  R.  Ercoline,  and  H.  L.  Hastings  (Ref  15)-  This 
program  was  used  to  find  distribution  functions  for  compar- 
ison to  the  analytic  solutions,  to  investigate  how  to  best 
fit  the  idealized  analytic  solutions  to  real  phenomena,  and 


to  help  verify  the  limits  of  applicability  of  the  analytic 
solutions . 


I 

i 

( 

I 
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II.  Solution  of  the  Time  Independent 


Boltzmann  Eg uation  for  Inelastic  Collisions 


Rationale  for 


The  published  literature  contains  examples  that  indicate 
that  F(€)  approximates,  at  high  enough  energy,  a Maxwellian 


distribution 


1 

exp^j^ 


(see  Fig  1)  ( Refs  3;  and  12).  If 


the  fluctuations  from  an  ideal  Maxwellian  are  slowly  varying 
functions  of  energy,  it  seems  plausible  to  represent  F(€)  by 
a near-Maxwelli.an  ~ exp  . Now  T is  no  longer  a 

true  temperature  characteristic  of  the  distribution  as  a whole, 
but  is  a slowly  varying  function  of  €.  Hopefully  this  repre- 
sentation will  result  in  a form  of  the  Boltzmann  equation  that 
can  be  readily  solved. 


Derivation  of  the  Equation  f or  T(g ) 

Boltzmann  Equation  for  F(€ ) . It  is  assumed  that  the 
electrons  are  in  the  presence  of  a background  gas  of  identical 
massive  particles  (which  may  be  atoms  or  molecules).  The 
particles  have  one  or  more  excited  states;  the  particles  are 
massive  enough  that  elastic  collisions  may  be  ignored.  The 
electron  density  is  assumed  to  be  low  enough  for  electron- 
electron  interactions  to  be  ignored.  It  is  also  assumed  that 
the  environment  (and,  as  a consequence,  F(€))  is  spatially 
homogeneous  and  isotropic  in  energy  space.  These  assumptions 
are  quite  reasonable  for  a typical  electric  discharge  laser. 
Ionization  will  not  be  considered. 

Under  these  conditions  equation  (1-2)  can  be  derived 
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f(€)  (€“ 


Figure  1.  Published  Distribution  Functions 
in  N2  (Ref  12:  1992) 


I 


from  Rockwood’s  formulation.  The  reverse  reaction  rate  con- 


i 


I 

I 


I 

I 


stants  are  related  to  the  forward  rates  by  the  principle  of 
micro-reversibility  (Ref  2:  431). 

Form  for  £ and  R (€ ) . It  would  be  extremely  difficult  to 
solve  equation  (1-2)  analytically  using  the  experimentally 
derived  rate  constants  and  an  electron-energy-dependent  €. 

Therefore,  some  simplifying  assumptions  will  be  made  about 

• • 

€ and  the  form  of  R(€).  First,  6 is  assumed  to  be  constant 
in  energy  space.  Using  an  expression  for  € derived  from 
Rockwood,  it  can  be  shown  that  this  assumption  is  equivalent 
to  assuming  that  the  collision  frequency  of  electrons  with 
the  background  gas  is  independent  of  energy  (Ref  13!  2348). 

Although  this  assumption  only  approximates  the  real  cases, 
it  does  simplify  the  solution  of  the  problem  considerably. 
Second,  the  following  form  for  R0(^)  is  assumed: 


R3(€)  = 


0,  € < Xg 

Rp.  € > 


where  Xq  is  the  threshold  energy  for  the  excitation  to  the 

i 

excited  state  3*  Since  R is  proportional  to  where  <r 

is  the  inelastic  collision  cross  section,  this  assumption  is 

_A 

equivalent  to  assuming  that  a is  proportional  to  € ® (see 
Fig  2).  Again,  this  assumption  only  approximates  real  cross- 
sections.  However,  it  simplifies  the  solutions  considerably, 
and  it  appears  from  the  numerical  calculations  discussed  in 
Chapter  III  that  the  solutions  are  insensitive  to  the  exact 
form  of  the  cross-sections. 
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1 


Figure  2.  Analytic  R and  a 


Derivation  of  Eg uation  for  T(g ) . For  convenience  in  the 
derivation  and  in  comparisons  of  analytic  results  to  numerical 
results,  a function  f(€)  was  defined  such  that 
f(€)  = €^F(€)  . Consistent  with  the  assumption  that 

F(e)  ~ e*  exp|j^]  for  high  energies,  an  assumed  form 

f(€)  = c |exp  ^ p was  used.  Finally,  from  published  nu- 

merical results,  it  was  observed  that  f(C)  drops  off  over  an 
energy  range  much  less  that  €,  and  does  so  fairly  smoothly, 
for  € » X , as  long  as  E/n^  is  less  than  approximately 
4 X 10~^^  Vcm^  (See  Fig  1)  (Refs  3?  and  12).  This  led  to 
two  more  assumptions,  which  are  crucial  to  the  derivation: 


12 


(which  implies 


1)  » 1 , for  € > Xp 


2)  -i-  ^ ^ for  € > X 

kT  ae  ^^2  ’ ■^a  • 


g»i ) 


(2-2) 


In  Appendix  A,  it  is  shown  that  these  assumptions  are  consis- 
tent with  the  solution  found. 

Using  these  definitions  and  assumptions,  a straightfor- 
ward but  tedious  derivation  leads  from  Eq  (1-2)  to 


3(kT)^  3k 


-66^  ^ 2gg^/aT\  ^ . 'Sp  r / g 

.2,j,3  ae  g ®^P^kTTr 

- Vb<^>1  • '=■='8  '2-3> 


Solution  for  f (€  ) , Two-Level  Case 

€ > X.  For  the  two-level  case  being  considered,  3 can 
represent  either  the  ground  state  or  the  excited  state. 
Since  €’  is  the  electron  energy  of  the  initial  state  of  a 
reaction  where  the  electron's  final  energy  is  €, 


e*  = 


€ - Xg  , for  3 = ground  state 
€ + Xg  , for  3 ” excited  state 


(2-4) 


From  the  principle  of  micro-reversibility,  it  can  be  shown 
that  Rp(€)  = Rp(€  + Xp)((€  + Xg)A)=  and 

Rp(e  - Xp)  = Rp(€)(€/(e  - Xg))=  (Ref  2:  431).  Using  n^  for 

the  ground  state  particle  density,  n^  for  the  excited  state 
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L 


particle  density,  and  the  assumed  values  for  results  in 


0 = 


dT  ^ 2e€^/aT\^  ^ / g 

3(^kT)2  ■ jk^'P  ®^p(ETm 


2€€ 


T^e'-^X))  ■ ■*■  ®^P(kT{€) 


kfT 


(6  + X)  \ /g  + X\^  o,,  ^ /€  + X\^ 

(2-5) 

where  the  "3"  subscripts  have  now  been  dropped  and  micro- 
reversibility has  been  used.  Using  the  assumptions  already 
made,  and  neglecting  super-elastic  collisions  (which  follows 

from  the  assumption  that  n « « ),  it  is  straightforward 

X o 

to  show  that,  for  € > X 


0 = 


2g€ 


3(kT)‘ 


466^  dT  _2iil/lT\^ 
3k2^3  3^2^^ lag/ 


g > X 

(2-6) 


This  equation  has  the  general  solution 


“ k (\^) 


Ve  + c' 


g > X 


(2-7) 


which  leads  to  the  general  solution  for  f(g)  of 


/6n  Rg 

f (g  ) = c exp  {-■'  ° 


g > X 


(2-8) 


where  c is  a constant  of  integration,  which  will  be  evaluated 
later.  For  convenience,  the  constant  factor  in  the  exponen- 
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tial  will  be  defined  as  follows: 


so  that 


f(€)  = cfexp 


^exi 


€ > X 


(2-9) 


(2-10) 


€ < X.  Since  f(e)  is  now  known  for  £ > X , the  coupled 
terms  in  the  collision  term  of  equation  (1-2)  can  be  evalu- 
ated.  If  the  previous  assumptions  are  used,  and  €®f(€)  is 
substituted  for  F(€),  equation  (1-2)  eventually  becomes 


i = /■('vF-nt)  exp|-^  (2-11) 

where  c is  the  constant  of  integration  in  Eq  (2-8). 

The  integral  in  equation  (2-11)  can  be  evaluated,  but 

the  resulting  expression  for  ^ cannot  be  evaluated  in  closed 

form.  Therefore,  the  integral  in  equation  (2-11)  will  be 

approximated  at  € * X and  6 « X . The  two  solutions  for 

f(€)  that  result  from  integrating  these  approximate  expres- 
d f 

sions  for  ^ will  be  matched  at  the  boundary  between  them. 

For  € * X , the  upper  limit  of  the  integral  in  equation 
(2-11)  is  set  equal  to  X.  Integrating  twice  leads  to 

f(e)  = c expl^l  -4 + c,  , q < € < X (2-12) 

^^^x  * e^kT^  ^ 

where  c^^  is  a constant  of  integration  and  q is  the  boundary 


between  the  two  solutions  for  € « X and  € » X 

Jx  •Jv  + X 

With  the  change  of  variable  u = ^ , equation 


(2-11)  may  be  rewritten  as 


M. 

36 


2€*  VX 


€ < X 


(2-13) 


For  e « X , the  upper  limit  may  be  approximated  by 

X /.  ■ € 

kT 


(l  + . If  the  integral  is  then  evaluated,  using 


X/kT  » 1 , the  result  is 

3C 


3€ 


-c  exp 


-X 

kT. 


2kT 


© ■ ®^P[2kT  ]}  ' ^ (2-14) 


This  cannot  be  readily  integrated.  However,  by  using  the 
McLaurin  series  for  the  exponential,  the  equation  can  be  re- 
duced to 


df 

-c  exp  , 

/X\' 

CD  1 

z 

-1 

d€ 

2kT  \e/ 

2kT„ 

X 

II 

L xj 

For  € = q 

(the  worst  case). 

and 

q = k 

n 


-n 


n 


!(n  -^T 


€ < q 
(2-15) 


discussed  in  Appendix  B,  the  first  three  terms  A^  in  the 
series  are 
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= 1 


Ag  = .0833 

A^  = .00833 


(2-16) 


It  should  be  obvious  that  retaining  the  first  two  terms  alone 
gives  excellent  results.  The  resulting  equation  may  be 
integrated  to  give 


f(6)  = 


c exp  (-X/kT^) 


^ < q 


3(kT^)^ 


(kT^)' 


(2-17) 


where  is  a constant  of  integration. 

F ( € ) 

Complete  Solution.  Since  ^ “^dg  represents  the  propor- 
tion of  the  total  electron  population  with  energies  € to 
€ + d€  , / ^ dg  should  equal  1.  This  leads  to  the 

In 

normalization  condition 


^f  (g)dg  _ 


(2-18) 


Utilizing  Eq  (2-17),  q = kT^  , continuity  of  f(g),  and 
defining  f(X)  = 1 , the  following  solution  is  found; 


f(€)  = exp 


(A\1 
1 - 

[x)  ■ 


€ > X 


(2-19) 


f(€)  = 1 ^ i^KD"  - l]  . kT^<  € < X 
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f(e)  = 1 + 


kT^ll2\kTy_/ 


1 + 


6 < kT 


^"2  \j_  + + 

^ [18  3kT 


(2-19) 


2kT  4(kT  4(kT 


By  defining  a dimensionless  variable  ^ ~ x ^ parameter 

X 

® = ^j-Ts—  , the  solution  can  be  written  in  the  following  form: 
■^x 


f(u)  = exp(a(l  - ^u))  , u > 1 


f(u)  = 1 + a(v/l/u  - 1)  , l/a  < u < 1 


(2-20) 


f (u)  = 1 + a I^n/^  ~ ^ - 2uM  , 0 < U < i 


N = X'*  JL  + “ + 2 ^ JL  + JL 
^ ^ 18  ^ 3 « «2  \,3 


It  is  important  to  note  that  the  form  of  the  solution  is 
determined  entirely  by  a . This  is  discussed  in  more  detail 
in  Chapter  III. 


Solution  for  f (€  ) , Three-Level  Case 

Nuw  a more  complicated  case  will  be  considered.  The 
background  gas  will  be  assumed  to  have  two  excited  levels, 
with  thresholds  and  X^,  and  rate  constants  and  R2  above 
threshold.  It  will  be  assumed  that  X^  < (X2  " X^^)  . This 

assumption  is  made  for  two  reasons.  First,  it  approximates 
fairly  well  a large  number  of  molecular  systems  that  have  one 
or  more  vibrational  levels  at  fairly  low  energies,  and  then 
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electronic  excitational  levels  at  higher  energies.  The 
approximation  is  especially  good  if  the  spacing  among  the 
vibrational  levels  and  among  the  electronic  levels  is  small 
compared  to  the  spacing  between  the  vibrational  levels  and 
the  electronic  levels.  Second,  the  assumption  makes  the  form 
of  the  solution  much  simpler.  It  should  be  noted  that  this 
assumption  is  not  strictly  necessary  for  the  solution  of  the 
problem.  However,  without  this  assumption  the  complexity 
of  the  problem  increases  greatly  without  concomitant  increases 
in  physical  understanding. 

€ > X.  For  both  € > X2  and  < € < X2  it  is 
straightforward  to  show  that  equation  (1-2)  leads  to  equa- 
tions formally  identical  to  the  two  level  case.  The  re- 
sulting solutions  are 


f(€)  = c_  exp 


-4 

L “xaJ  ' 

r ^/X,€  -I 

f(e)  = exp[^  J , 


€ > X, 


Xi  < e < X2 


(2-21) 


where 


yxT 


kT. 


x2 


kT. 


6n  R, 
o 1 


xl 


It  is  still  necessary  to  assume 

n n 

» 1 , and  — ^ » 1 . 


kT, 


xl 


» 1 , 


kT 


» 1 


x2 


n 


xl 


n 


xl 
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€ < It  is  also  straightforward  to  show  that  equa- 

tion (1-2)  leads  to  an  equation  for  af/ae  that  has  a sum  of 
terms,  each  of  which  is  formally  identical  within  a multipli- 
cative constant  to  the  term  in  Eq  (2-].l).  The  resulting  so- 
lutions are 


f(€)  = -Y  exp  I 


xl-* 


i (rj  * rJ 


h 

^ ^ 


q < € < X. 


f(€)  = exp  ^-—3  - , ,2  } 


!2  expr-:^^! 


X| 

3(kT^j)3 


^4  ^ 


€ < q 


(2-22) 


Complete  Three-Level  Solution.  By  defining  f(X2)  = 1 1 
using  continuity  of  f(€)  and  Eq  (2-18),  and  choosing 
q = kT^^  (this  choice  is  discussed  in  Appendix  B),  the  fol- 
lowing expressions  for  f(€)  and  N are  developed  for  the  three- 
level  model. 


f(€)  = exp  -X 


f(€)  = exp 


. X2  < € 


. X,  < e < x^ 
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[Rl  + Rg)  kT^2 


) ■ 


kT  , < € < X, 
xl  1 


■ -[^1 4?  -^) 


^7  I b-  /F\ 

* (Hi  * Rj)  (aSkT^  - t) 

/iz/iLt . A 


* 2(1*1  " «2>U(M^3)2  ■ 12(M^p3  * “x2<“xl>‘'' 


^^(Xi)'^}  ■ ^ < 


(2-23) 


« (r-,  1 2(kT  )3  ) 2(kT,  )2Vx- 

N = 2{>/}C(kT  ^ — f 7-2^^ 2 ^ xl 2 


2'  x2' 


2(kT  )^)  - X^1  ( _ 

* ^ Jl_  ■*, 

x4  6>*Txl  36  ) 5(R~r-F 


4(kT^^)2 


R^T  kT^2  \ h - 


12(kT^2^^  ^^^x2 
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Extension  of  the  Solution 


More  Levels  ^ ^2  ~ Xj^) . Although  the  algebra  would 
become  quite  involved,  the  extension  of  this  model  to  four 
or  more  levels  would  be  quite  straightforward,  as  long  as 
the  appropriate  assumptions  about  and  ^ can  still  be 
made.  As  long  as  € is  greater  than  some  threshold  energy 
f(€)  will  have  the  form  of  equation  (2-20),  with  R re- 
placed by  a sum  of  over  all  states  3 with  threshold  ener- 
gies less  than  or  equal  to  X^^.  For  3 less  than  any  threshold 
energy,  the  differential  equation  for  3f/a€  can  be  reduced 
to  the  form  of  a sum  of  terms,  each  of  which  is  formally 
identical  within  a multiplicative  constant  to  the  terms  in 
equation  (2-11).  This  results  in  forms  similar  to  equation 
(2-23),  but  even  more  complex. 

X^  > (X^  - Xj^).  (See  Figure  3). 


^ f- 


X2-X1  X2 


a.  X^^  < (Xg  - X^) 


r "i-r—  "2  — r ^3  ^ 


1 ^2 

€ 


b.  X^  > (Xg  - Xj^) 


Figure  3*  Comparison  of  threshold  energy  relationships 


The  extension  in  this  case  would  be  fairly  straight- 


forward. The  distribution  function  at  some  energy  € is 
determined  by  the  rate  electrons  are  lost  at  that  energy  due 
to  inelastic  collisions,  if  those  collisions  are  possible 
at  €.  This  is  the  process  applicable  for  f^  and  f^  in  Fig 
3a,  and  f^  and  f^^^  in  Fig  3b.  If  € is  less  than  all  the 
threshold  energies,  then  inelastic  collisions  cannot  take 
place.  In  this  case  the  distribution  at  € is  determined  by 
the  rate  at  which  electrons  are  gained  from  inelastic  colli- 
sions at  energies  f , which  is  proportional  to 

f(€  + Xp)  . This  is  the  process  applicable  for  f^^  in  Fig  3a, 
and  fj^  and  in  Fig  3b.  In  Fig  3a,  f(€  + X^^)  has  the  same 
form  for  any  energy  € in  the  range  0 < t < Xj^  . As  a re- 
sult, one  differential  equation  is  valid  for  the  entire 
range.  In  Fig  3h,  f(€  + X^^)  changes  form  as  € goes  from 
0 < 6 < (X^  ~ "to  X2  - Xj^  < € < X^  (or  as  € goes  from 

fj^  to  ^2^’  ^ result,  two  differential  equations  are  neces- 

sary for  the  range  0 < € < X^^.  This  does  not  cause  any  funda- 
mental difficulties,  but  does  lead  to  a much  more  complicated 
expression  for  f(€). 


Ill  Analysis  and  Conclusion 


General  Considerations 

The  analytic  two-level  solution  has  two  parameters;  X 

and  kT^.  However,  it  is  possible  to  express  f(€)  in  terms  of 

X € 

and  the  dimensionless  parameter  • As  a result, 

it  is  more  appropriate  to  characterize  the  various  solutions 
by  the  corresponding  values  of  X/kT^:  a].l  the  solutions  pre- 
sented here  are  so  identified. 

Much  the  same  situation  exists  for  the  three-level  case, 
except  that  the  parameters  are  now  X^ , X^. 

Rl  . 

P — - . With  this  many  parameters,  it  would  be  impracti- 

cal to  try  to  represent  all  possible  combinations,  and  there 
has  been  no  attempt  to  do  so  here. 

All  the  cross-sections  used  for  comparison  of  the  an- 
alytic solution  to  realistic  gases  are  from  Kiefer  (Ref  6). 
Since  each  cross-section  will  be  identified  at  the  time  it  is 
introduced,  there  will  be  no  further  citations  of  the  cross- 
section  source.  All  numerical  calculations  were  done  with 
the  program  NGB,  as  discussed  in  Appendix  C. 

Criteria  for  Evaluation 

A discussion  of  the  criteria  used  to  evaluate  the  ana- 
lytic solutions  developed  in  this  thesis  would  be  in  order 
at  this  time.  In  general,  analytic  models  cannot  describe 
physical  situations  as  complicated  as  the  one  involved  here 
as  accurately  as  numerical  solutions.  Therefore,  it  would 
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f 


I 

f 

( 

be  inappropriate  to  use  exact  agreement  with  the  numerical 
I solution  as  a general  criterion  for  evaluating  the  numerical 

f solution,  although  such  agreement  certainly  adds  support  to 

; the  validity  of  both  the  numerical  and  analytic  solutions. 

t 

I More  appropriate  criteria  would  be  those  that  measure  the 

[ insight  the  analytic  model  gives  into  the  physical  situation. 


i 

i 

i 


For  instance,  it  is  more  important  to  be  able  to  demonstrate 
the  dependence  of  f on  € than  to  produce  hard  numbers.  Based 


on  this  principle,  the  following  criteria  were  used  to  eval- 
uate the  validity  and  utility  of  the  analytic  solutions : 

1.  Does  the  analytic  solution  dependence  on  electron 
energy  agree  with  that  of  the  numerical  solution, 
within  the  limits  of  the  numerical  solution? 

2.  Can  the  analytic  solution  be  used  to  predict  the 
scaling  of  physical  parameters  in  a meaningful 
fashion? 

Obviously,  the  first  condition  is  met  for  those  cases 
where  the  two  solutions  are  essentially  indistinguishable. 


but  even  if , there  are  differences  between  the  two  solutions. 


this  condition  can  still  be  met  if  the  two  solutions  have 
similar  slopes  at  the  same  energy  (See  Fig  9)- 

As  long  as  the  first  condition  is  met,  the  analytic 
solution  can  be  considered  valid.  It  can  be  considered  use- 
ful if  the  second  condition  is  met.  The  success  of  this 
model  in  meeting  the  second  criterion  will  be  discussed  later 
in  this  chapter. 

Of  course,  it  is  possible  that  an  analytic  model  will 
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demonstrate  other  results  in  excess  of  the  criteria  listed 
above.  Such  serendipitous  results  may  increase  the  utility 
of  such  a model,  but  the  lack  of  them  does  not  invalidate 
it. 

Theoretical  Cross-Sections 

Two -Level  Model . Figures  4, 5, 6, and  7 show  the  analytic 
and  numerical  solutions  for  values  of  X/kT^  of  10,  20,  30, 
and  100.  The  cross-sections  used  were  calculated  to  produce 
R's  constant  above  threshold.  The  background  gas  contained 
no  excited  atoms  and  the  numerical  calculations  included  only 
the  effects  of  heating  and  inelastic  collisions. 

Several  points  should  be  made.  First,  the  apparent  dis- 
crepancy at  the  extreme  high  energy  end  of  the  energy  axis 
is  purely  a result  of  NGB.  NGB  always  conserves  the  total 
electron  density.  As  a result,  electrons  being  driven  up  the 
energy  axis  by  heating  must  eventually  stop  at  the  upper 
boundary.  This  causes  the  increase  in  f(€)  at  the  upper  end 
of  the  axis.  (In  some  cases  the  entire  energy  axis  used  in 
NGB  is  not  depicted,  in  which  case  this  discrepancy  is  ab- 
sent). 

Second,  there  is  a discrepancy  at  the  low-energy  end  of 
the  energy  axis.  This  is  also  a result  of  the  finite 
differencing  used  in  NGB.  The  behavior  of  f(e)  for  € * 0 
is  strongly  influenced  by  f(€  + X).  The  distribution  func- 
tion drops  very  rapidly  just  above  X.  If  the  resolution  of 
NGB  is  larger  than  the  interval  over  which  f(€)  drops,  it 
will  be  unable  to  accurately  describe  the  behavior  at  € ® X, 
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Figure  4.  Ideal  Two  Level  System, 
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and  also  at  € = 0 as  a result  of  the  strong  coupling.  In- 
creasing the  resolution  in  NGB  would  improve  the  agreement 
between  the  analytic  and  niamerical  solutions. 

Figure  8 shows  two  numerical  solutions  and  the  analytic 
solution  for  X/kT  = 15  ■ Notice  that  as  the  energy  reso- 

X, 

lution  increases,  the  discrepancy  at  the  low  end  of  the  en- 
ergy axis  decreases.  Note  also  that  as  the  upper  end  of  the 
energy  axis  moves  dov/n,  the  discrepancy  at  that  end  moves 
down  also.  This  clearly  demonstrates  that  the  discrepancy 
at  that  end  is  a result  of  NGB  and  does  not  reflect  physical 
reality . 

Finally,  it  is  apparent  that  the  numerical  and  analytic 

y 

solutions  agree  very  well,  especially  foi'  values  of 

■^x 

greater  than  ~ 15 

Three-Level  Model . Figure  9 shows  analytic  an  numer- 
ical solutions  of  an  idealized  three  level  system.  The  nu- 
merical and  analytic  solutions  do  not  agree  as  well  for  this 
case  as  for  the  various  two-level  cases,  especially  for 
€ * . Howei^er,  note  that  the  scaling  for  e > X^  is 

still  very  good,  as  evidenced  by  the  fact  that  the  two  curves 
have  very  similar  slopes.  In  addition,  a comparison  of  the 
average  pumping  rate  for  the  second  level,  given  by 


J'R^{^)F{€)d€ 


(3-1) 


results  in  a value  of  = 6.9  x 10”^®  cm^/sec  for  the  an- 

__  1 O 

alytic  case,  versus  3*6  x 10  for  the  num.erical  (Equation 
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(3-6)  gives  niore  explicitly).  As  described  in  Nielsen’s 

unpublished  paper,  the  pumping  rates  of  the  higher  levels  can 
be  very  sensitive  to  the  width  of  the  energy  bins.  Since  f(€) 
drops  off  rapidly  within  a range  ~ kT^  above  threshold, 
most  of  the  electrons  that  lose  energy  do  so  with  a range 
~ kT^  above  that  threshold.  These  electrons  are  spread 
throughout  the  first  energy  bin,  if  that  bin  is  wider  than 
kT^.  This  error  in  the  electron  density  in  energy  space  can 
cause  errors  throughout  the  energy  axis,  due  to  the  strong 
coupling  of  the  Boltzmann  equation.  In  the  example  shown  in 
Figure  d6  was  .1  eV  , while  kT^j^  was  .129  eV,  and  hT^2' 
,176  eV.  With  d€  so  close  to  kT^^,  it  is  difficult  to  say 
which  of  the  two  average  rates  is  more  accurate:  computer 
limitations  made  it  impractical  to  use  higher  resolution  on 
the  energy  axis.  In  any  case,  the  agreement  is  certainly 
close  enough  to  be  useful  in  predicting  realistic  reaction 
rates,  as  will  be  seen  later. 

Note  that,  unless  6 is  very  large,  » <R2>  . 

This  means  that,  as  for  the  two-level  case,  conservation  of 
energy  implies 

^^1^%  " X (3-2) 

independent  of  the  exact  form  of  the  distribution  function, 
as  long  as  <Ri>  » ^R2^  • For  this  reason  is  not  a 

useful  test  for  the  two-level  model,  nor  is  ^Rj^  for  the 
three  level  model. 

I 
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k Realistic  Cross-Sections 


Two -Level  Model . Figures  10,  11,  and  12  show  compari- 
sons of  the  analytic  and  numerical  solutions  using  realistic 

cross-sections.  Figure  10  is  based  on  the  cross-sections  for  = 

i 

the  lowest  electronic  state  for  N2.  Figure  11  is  based  on  the  | 

second  electronic  state  for  N^.  Figure  12  is  based  on  the  ' 

second  electronic  state  for  Ar.  Several  different  methods 

of  matching  the  analytic  cross-sections  and  thresholds  to 

the  realistic  ones  were  used.  The  method  finally  used  was  to 

choose  an  analytic  R of  one-half  the  maximum  reaction  rate  for 

the  real  cross-section,  where  the  reaction  rate  is  given  by 

'*<«>■■  (I;)” <3-3) 

» 

where  m is  the  electron  mass 

a is  the  inelastic  collision  cross-section  (Ref  13: 

23^9). 

j 

I 

The  threshold  was  chosen  equal  to  the  energy  lost  by  the  • ; 

electron  upon  collision.  This  method  of  matching  was  used 
primarily  because  it  worked.  This  point  will  be  discussed 
in  more  detail  later. 

From  Figures  10,  11,  and  12  it  seems  that  the  analytic 
and  realistic  solutions  match  fairly  well,  except  at  the 
high  energy  end  of  Figures  10  and  12.  This  discrepancy  is 
caused  by  the  exact  behavior  of  the  realistic  cross-sections. 

The  analytic  model  assumes  that  there  is  a finite  cross- 
section  for  all  energies  above  the  threshold.  In  this  case, 
however,  the  real  cross-se  tions  go  to  zero.  As  a result 
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(since  elastic  collisions  are  excluded),  there  is  no  mech- 
anism to  allow  the  electrons  to  lose  energy,  and  the  dis- 
tribution approaches  an  infinite-temperature  Maxwellian 
distribution. 

Three-Level  Model . Figures  13 > l^i  and  15  compare  the 

numerical  and  analytic  results  for  three  different  cases. 

Figure  13  is  based  on  the  first  and  fourth  electronic  levels 

on  Ng*  Figure  14  is  based  on  the  second  and  fifth  electronic 

levels  of  Ng.  Figure  15  is  based  on  the  first  and  fifth 

levels  on  N^.  In  all  cases  the  analytic  reaction  rates 

were  chosen  to  be  one  half  the  maximum  realistic  rates j the 

thresholds  were  chosen  equal  to  the  energy  loss  for  each 

transition.  The  agreement  in  Figures  13  and  14. is  fairly 

close;  however,  the  agreement  in  Figure  15  is  poor.  The 

analytic  and  numerical  pumping  rates  for  the  second  excited 

level  for  Figure  13  are  2.78  x 10  cm^/sec  and 

6.42  X 10  cmvsec  , respectively.  For  Figure  14,  the 

values  are  4.42  x 10~^^  cm^/sec  and  4.53  x 10"^^  cm^/sec  , 

respectively.  For  Figure  15,  the  values  are  8.49  x 10~^^ 

3 “113 

cmVsec  and  1.14  x 10~  cm^/sec  , respectively.  Again, 
the  results  in  Figures  13  and  14  agree  well:  those  in  Figure  | 

15  do  not.  The  discrepancy  is  caused  by  an  effect  similar  ii 

to  that  just  discussed  for  Figures  10  and  12.  For  the  data 
depicted  in  Figure  15,  there  is  a region  between  and  X2 
where  all  the  cross-sections  are  zero.  As  a result,  any 

electron  that  reaches  this  region  cannot  lose  energy  by  ' 

exciting  a background  atom.  Instead,  it  must  keep  increas- 
ing in  energy  until  it  reaches  X^.  This  causes  f(6)  to  be 
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Figure  13-  Three  Level  System, 
El  and  E4  Levels  of 


much  higher  in  this  region  than  the  analytic  model  predicts. 
In  the  conclusions,  a method  to  correct  this  discrepancy  will 
be  discussed.  On  the  other  hand,  if  the  two  cross-section 
profiles  overlap,  the  process  described  above  does  not  occur, 
and  the  result  agrees  with  the  analytic  model.  Although  this 
effect  limits  the  number  of  real  cases  which  can  be  described 
by  this  model,  it  does  not  significantly  reduce  the  effec- 
tiveness of  the  model  in  meeting  the  criteria  discussed  ear- 
lier, especially  since  it  is  possible  to  predict  a priori 
when  this  effect  will  occur. 

Modeling  a Multi-Level  System . Figure  16  shows  a com- 
parison between  a five-level  system  and  a three-level  an- 
alytic model.  The  five-level  system  is  intended  to  repre- 
sent a system  with  characteristics  similar  to  a real  system. 
The  system  uses  cross-sections  for  the  sixth  and  seventh 
vibrational  and  the  first  and  second  electronic  levels  of 
CO.  Program  limitations  in  NGB  made  it  impractical  to  use 
all  seven  vibrational  levels.  These  limitations  could  be 
easily  removed,  if  necessary.  However,  in  this  particular 
case,  such  modification  would  not  be  worth  the  effort  in- 
volved. For  the  analytic  model  used  one  half  the  av- 
erage of  the  maximum  rates  for  the  V6  and  V?  levels.  For 
R2.  the  analytic  model  used  one  half  the  average  of  the 
maximum  rates  for  the  El  and  E2  levels.  For  X^,  the  analytic 
model  uses  the  average  of  the  energy  losses  for  the  V6  and 
V7  transitions.  For  the  analytic  model  uses  the  average 
of  the  energy  losses  for  the  El  and  E2  levels.  The  numer- 
ical results  for  the  pumping  rates  for  El  and  E2  are 
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= 2.07  X 10“^®  cmVsec  and  ^^^2)  ~ ^ lO"^^  cmVsec. 

The  analytic  rates  for  El  and  E2  are  given  by  Equation  (3-6) 
with  and  replaced  by  and  as 

appropriate,  and  with  X^,  and  Rj^  as  described  above.  The 
results  are  = 2.6  x 10~^®  cm^/sec  and 

^Re2^  = 2.6  X 10  ^ cm^/sec  . The  poor  agreement  for  the  E2 
level  is  not  surprising:  the  analytic  model  assumes  the  level 
involved  is  well-separated  from  any  lower  thresholds,  a condi- 
tion which  certainly  does  not  hold.  Since  f(e)  drops  rapidly 
above  threshold,  the  presence  of  a higher  threshold  is  not 
as  significant.  The  close  agreement  for  the  El  level  and 
for  the  distribution  function  itself  indicates  that  the  an- 
alytic solution  would  be  useful  for  modeling  cases  similar  to 
this  one. 

Superelastic  Collisions . Figure  1?  shows  the  result  of 

a numerical  and  analytic  comparison  used  to  investigate  the 

validity  of  neglecting  superelastic  collisions.  The  analytic 

model  used  a value  of  X/kT^  of  25.  The  numerical  calculations 

used  idealized  reaction  rates  constant  above  X = 5 eV 

According  to  the  expression  -developed  in  Appendix  A (Eq 

(A-11))  superelastic  collisions  should  be  negligible  as  long 
/ 4 

as  n^/n^  » 3*1  x 10  . The  numerical  results  shown  had 

values  of  n /n  of  lo'^,  10^,  and  10^  . Values  of  n /n 

OX  ox 

7 

greater  than  10  produced  results  indistinguishible  for  the 
result  with  n^/n^  = 10'  . As  can  be  seen,  the  result  starts 

to  change  for  n /n  = 10^  . For  n /n  = 10^  , the  numer- 

OX  OX 

ical  solution  has  values  much  higher  than  the  analytic  at 
high  energies,  and  lower  values  at  low  energies.  This  is 
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exactly  the  expected  result.  Superelastic  collisions  act  to 
transfer  electrons  from  low  energies  to  higher  energies;  if 
this  transfer  involves  significant  numbers  of  electrons,  the 
effect  is  exactly  that  produced  by  NGB.  Note  that  the  thres- 
hold for  this  effect  to  become  significant  agrees  very  well 
with  that  predicted  by  Eq  (A-11). 


Conclusions 

For  a system  composed  of  an  electron  gas  in  the  presence 
of  background  particles  capable  of  being  excited  internally 
by  the  electrons,  the  form  of  the  electron  distribution 
function  depends  strongly  on  the  parameters  of  the  form 


(3-4) 


where  n^  is  the  background  gas  particle  density 
R is  the  average  reaction  rate 
X is  the  excitation  threshold 

€ is  the  average  electron  heating  rate  (per  electron) 

X ~ 

as  long  as  > 10  and  elastic  and  superelastic  collisions 

can  be  ignored. 

For  heating  produced  by  a uniform  non-time-varying  elec- 
tric field,  this  result  reduces  to 


^ - E 

X ~ n V^RXmv 
o ^ 


(3-5) 


where  E is  the  magnitude  of  the  electric  field 
e is  the  electron  charge 
m is  the  electron  mass 


* 
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V is  the  average  electron-background  particle  colli- 
sion frequency 


For  typical  systems, 


» 1 implies  — « 10“^^ 


For  laser  media  under  such  circumstances,  if  the  excited 
particle  density  is  low  enough  that  superelastic  collisions 
can  be  ignored,  and  if  the  background  gas  can  be  approxi- 
mated as  having  two  excited  levels,  the  electron  pumping  to 
the  second  level  is  given  by 


<^z> 


,2  R„X„e  n^ 
2 2 o 


^ - ^il  A 

n^  V6R,mV 
o i 


^^l^^l  " ^2- 


(3-6) 


" ^2 


y‘/«.  y 
^2  ^1 


where  use  has  been  made  of  » 1 and  ^2  ^ ^1 

The  importance  of  the  parameter  x/kT^  and  the  manner  in 
which  it  affects  the  distribution  function  indicates  that, 
for  a numerical  solution  to  be  valid,  it  must  have  a resol- 
ution much  better  than  kT  ; for  example,  kl  /d€  > 10 

It  is  difficult  to  compare  this  conclusion  with  pub- 
lished results,  as  very  few  authors  give  the  resolution  used 
in  their  numerical  results.  It  is  interesting  to  compare 
with  the  results  that  are  available.  Rockwood  decreased  the 
energy  difference  in  his  method  until  further  reductions  had 
no  effect  on  the  result.  The  energy  difference  d€  chosen  in 
this  manner  was  0.01  eV.  Although  it  is  not  possible  to  cal- 
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culate  kT^  exactly  from  his  data,  it  would  be  close  to  0.005  eV, 
which  certainly  is  of  the  same  order  as  d€ . (Ref  13)  Nielsen 
et  alia . in  their  work  on  laser  induced  breakdown,  used  the 
photon  energy  of  the  laser  photon  as  d€ . In  their  CO^  laser 
breakdown  studies,  d€  = .1  eV  , and  kT  =*  .05  eV  (Ref  10). 
Again  the  agreement  is  fairly  good.  However,  in  the  ruby 
laser  breakdown  studies,  d€  = 1.7  eV  , and  kT^  .01  eV 
These  results  seem  to  invalidate  either  their  results  or  the 
analytic  model.  However,  note  that  the  analytic  model,  by 
its  choice  of  form  of  the  Boltzmann  equation,  assumes  a 
continuous  heating  process.  This  is  completely  valid  for  a 
D/C  field,  as  was  used  in  Rockwoods  research,  and  is  a very 
good  approximation  for  the  CO2  laser  in  helium,  where  the 
absorbed  photons  have  energies  of  0.1  eV  and  the  first  thres- 
hold is  ~ 20  eV  . But  for  a ruby  laser  in  deuterium,  where 
the  photon  energy  is  ^ 1.?  eV  , and  the  first  threshold  is 
still  ~ 20  eV  , it  is  difficult  to  justify  the  continuous 
heating  model . 

For  physical  situations  where  the  assumptions  the  analy- 
tic model  is  based  on  are  accurate,  the  analytic  model  des- 
cribes physical  reality  fairly  well.  In  particular,  it  is 

• Y 

necessary  that  € be  low  enough  for  the  assumption  » 1 ; 

that  be  small  enough  that  superelastic  collisions  can  be 
ignored;  that  n^  be  low  enough  that  electron-electron  inter- 
actions can  be  ignored;  and  that  the  heating  effect  be  con- 
tinuous in  energy  space.  Based  on  agreement  between  the  an- 
alytic and  numerical  solutions  discussed  previously,  the 
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final  solution  is  apparently  insensitive  to  the  exact  form 
of  R(€),  as  long  as  R(e)  / 0 for  € > X . As  a result,  it 
is  not  necessary  for  R(€)  to  closely  approximate  a constant 
value. 

The  analytic  two-level  solution  models  realistic  dist- 
ribution functions  based  on  the  two-level  systems  well,  espec- 
ially if  R(€)  is  non-zero  for  € > X . The  modeling  is  most 
accurate  for  kT  < € < X . Accuracy  of  the  model  for  € ^ 0 
is  difficult  to  determine  without  numerical  solutions  of  great- 
er resolution  than  are  available.  The  limiting  factor  in  such 
solutions  is  a combination  of  finite  computer  resources  and 
lack  of  resolution  in  experimental  cross-sections.  The  most 
accurate  method  of  approximating  R(€)  would  seem  to  be  to  use 
some  form  of  average  R similar  to 


where  "the  points  where  R(€)  is  0.  However,  in 

practice  this  method  is  not  clearly  adequate,  for  several 

^ reasons.  For  one  thing,  the  analytic  R extends  an  infinite 

! 

i distance  up  the  energy  axis;  this  makes  it  difficult  to  com- 

pare with  a realistic  rate  that  is  of  finite  extent  in  energy 
space.  Also,  any  such  average,  assuming  one  could  be  proper- 
ly picked,  would  generally  be  tedious  to  evaluate,  especially 
to  produce  what  is  at  best  an  estimate.  It  appears  that 
reasonably  good  results  can  be  acheived  by  using  R equal  to 
one  half  the  maximum  value  of  R(€)  for  the  particular  tran- 
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sition  involved,  and  by  choosing  X equal  to  the  energy  loss 
of  the  transition. 

The  analytic  three-level  model  can  be  used  to  model  real- 
istic three-level  systems  only  if  the  cross-sections  for  the 
two  excited  levels  overlap  and  the  condition  < (Xg  " 
is  met.  The  resulting  solution  can  be  used  to  predict  pump- 
ing rates  for  the  highest  level  within  an  order  of 
magnitude  or  less.  Good  results  can  apparently  be  obtained 
by  choosing  R21  X^,  and  Xg  the  same  as  in  the  two-level 
case . 

Under  some  circumstances,  the  three-level  system  can  be 
used  to  model  multi-level  systems.  The  thresholds  must  fall 
into  two  groups,  such  that,  if  Xj^  is  a representative  thres- 
hold for  the  lower  group,  and  X2  a representative  threshold 
for  the  higher,  then  X^  < (X2  " • There  must  be  no  in- 

tervals above  the  lowest  threshold  for  which  all  the  cross- 
sections  are  zero.  Good  results  can  apparently  be  obtained 
by  letting  the  reaction  rate  for  each  group  be  the  average  of 
one  half  the  maximum  reaction  rate  for  each  level  in  the 
group,  and  by  letting  the  threshold  for  each  group  be  the 
average  of  the  thresholds  for  each  group.  The  ptimping  rate 
for  the  lowest  transition  in  the  higher  group  can  apparently 
be  estimated  within  an  order  of  magnitude  or  less  by  an  ex- 
pression analogous  to  Eq  (3-6). 

It  would  be,  in  principle,  quite  straightforward  to 
extend  this  solution  to  more  levels,  or  to  levels  of  arbi- 
trary spacing,  and  such  extensions  would  allow  modeling  of 
cases  that  cannot  be  modeled  by  the  present  case.  However, 
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the  algebra  involved  in  such  extension  is  quite  involved. 

It  should  be  possible  to  extend  the  model  to  cases 
where  R(€)  above  threshold  is  zero  by  explicitly  including 
such  regions  as  separate  areas  for  which  a separate  solution 
must  be  found.  For  example,  if  R(€)  = 0 for  ^ ^2  ’ 

the  Boltzmann  equation  could  first  be  solved  for  € > 6^  ! 
that  solution  could  be  used  to  find  f(€)  for  ^ ^ ^ ^ ^ 2 ’ 
and  so  forth. 
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Appendix  A 


Validity  of  Assumptions 


This  assumption  may  be  validated  by  using  the  form  for 
T(€)  developed  in  Chapter  2.  It  will  be  evaluated  at  € = X , 
since  this  is  the  worst  case  condition.  From  such  an  evalu- 
ation the  following  expression  may  be  derived: 


€ « 6n  RX 
o 


(A-1) 


For  typical  values  of  n ~ 10^^  /cm^  , R ~ 10  cm^/sec  , 


and  X = 4 eV  (Ref  6),  this  gives 


€ « 2.4  X 10^^  eV/sec 


(A-2) 


This  value  may  easily  be  compared  with  published  results, 
From  Brown,  it  is  seen  that 


6 = 


2^2 
e E V 


2m 


(A-3) 


where 


e is  the  electron  charge 
m is  the  electron  mass 

V is  the  average  collision  frequency  of  the  electrons 
with  the  background  gas 
b) is  the  frequency  of  the  electric  field 
E is  the  amplitude  of  the  electric  field  (Ref  1:  l66) 
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For  a DC  field,  u)  = 0 . Also,  since  the  average  field 

strength  is  ^ times  the  amplitude  of  a sinusoidal  field,  but 

2 2 

is  equal  to  the  amplitude  of  a DC  field,  E /2  becomes  E . 

This  leads  to 


€ 


2 2 E 
n e — 
0 n 

g_ 

mi/ 


(A-4) 


1 2 

Using  the  same  typical  values  as  before,  and  v/  = 10  /sec 
(Ref  13:  2349)  gives 


^ « 3.7  X 10“^^  V-cm^  (A-5) 

This  is  just  above  typical  lower  values  of  E/n^  (Ref  12: 
1992-1993)-  This  result  is  a significant  limitation  on  the 
range  of  complete  validity  of  the  analytic  model.  Note,  how- 
ever, that  this  assumption  was  used  primarily  to  justify  neg- 
lecting the  contributions  at  some  energy  € (above  threshold) 
due  to  inelastic  collisions  at  an  energy  € > X.  If  the  as- 
sumption is  no  longer  valid,  then  those  contributions  can  no 
longer  be  neglected.  Even  so,  the  conclusion  about  the  re- 
solution required  in  numerical  solutions  are  still  valid. 

2 2 2 

For  fields  at  optical  frequencies,  u + v ^ u)  , and 
Eq  (A-3)  becomes 


. _ e E V 

o 2 

2mu) 


(A-6) 


Recall  that  the  average  Poynting  vector  amplitude,  S,  is 
given  by  (Ref  8:  465) 


Substituting  Eq  (A-7)  and  Eq  (A-6)  into  Eq  (A-1)  results  in 
6n  RXu>^c€m 

S « — (A-8) 

e V 

Using  the  same  values  as  bei'ore  and  u>  appropriate  for  a CO^ 
laser  results  in 

S « 1.2  X 10^  watxs/cm^  (A-9) 

This  is  near  or  even  above  reported  thresholds  for  breakdown 
in  clear  air  (Refs  10;  11).  It  is  certainly  high  enough  to 
allow  investigation  of  gas  laser  systems. 

_1_  dT  . - dfr 

Evaluation  of  this  expression  for  the  T(e)  developed  in 
Chapter  II  leads  to  a worst  case  of 

€ « 24n  RX  (A-10) 

o 

This  is  obviously  valid  if  Eq  (A-l)  is  valid. 

n « n 
X o 

This  assumption  is  used  on  several  occasions  to  justify 

discarding  terms  ~ n^f(€)R(€)  . The  worst  case  is  for 

C = 2X  . In  this  case  the  assumption  is 

n f(€  - X)  « n f(€)  . For  6 > 2X  , this  becomes 
X o 


5? 


elastic  collisions  could  be  ignored  in  mercury  vapor  for 
E/n^  greater  than  ~ 2 x 10  volts-cm  (Ref  13:  2353)- 

The  elastic  collision  rate  is  proportional  to  the  ratio  of 
the  electron  mass  to  the  background  particle  mass.  Mercury 
is  a very  heavy  particle,  with  a mass  some  4.5  times  heavier 
than  COgi  for  instance.  It  would  be  expected  that  Rockv/ood's 
case  would  represent  an  extreme.  Comparison  of  the  heating 
term  to  the  elastic  collision  term  in  the  Boltzmann  equation 
leads  to  the  requirement  that  X/kT^  < 100  for  elastic  colli- 
sions to  be  ignored,  which  is  not  inconsistent  with  the  re- 
sults above. 

The  best  approach  to  this  problem  would  be  to  repeat 
the  derivation  performed  in  this  thesis,  but  with  elastic 


collisions  included.  Unfortunately,  that  is  beyond  both  the 

58 


Appendix  B 
Choice  of  q 


General  Considerations 

Physical . Consider  that  f(€),  for  € < X , depends  on 
f(e  + X).  The  solution  for  € ~ 0 is  affected  most  strongly 
by  f(€)  for  € just  above  X (see  Fig  (B-l)).  Since  kT^  is 
characteristic  of  the  energy  range  for  which  there  is  a signi 
ficant  change  in  f(?),  for  6 > X , it  is  not  unreasonable 
to  expect  that  the  behavior  of  f(€),  for  € ~ 0 , might  also 

change  significantly  over  a range  ~ kT^ 


Figure  B-l.  Coir.parison  of  f(€) 


€ » 0 and  C « X 


^ Mathematical  There  are  several  approaches  that  could 

I 

be  used  to  choose  q.  One  would  be  to  choose  q so  that  N is 
insensitive  to  q:  that  is,  choose  q so  that  dN(q)/dq  = 0 
This  has  some  physical  significance  in  that  q is  strictly  a 
result  of  the  method  used  to  solve  the  Boltzmann  equation: 
its  effect  on  the  final  solution  should  be  as  little  as  pos- 
I sible. 

I Another  method  would  be  to  choose  q so  that  df/3€  is 

continuous  at  q.  This  approach  is  mathematically  satisfying, 
but  has  little  physical  meaning,  expecially  when  the  equa- 
; tion  for  df/3€  already  had  the  discontinuous  function  R(€). 

1 Regardless  of  the  approach  used,  q must  be  chosen  so  that 

[ 

r the  approximations  used  in  solving  the  equation  for  € < X 

: are  valid. 

( 

! 

Two-Level  Model 

dN/dq  = 0.  First  the  condition  dN/dq  = 0 was  used  to 
evaluate  q.  This  led  to  a complex  value  for  q,  which  was  un- 
acceptable. Since  it  was  impossible  to  find  a real  value  for 
q for  which  there  was  no  dependence  of  N on  q,  an  attempt 

was  made  to  find  a value  of  q which  minimized  the  dependence 

2 2 

of  N on  q by  using  the  condition  d N/dq  = 0 . This  led 

to  q = kT^ 

Continuity  of  df/d€  at  q.  In  order  to  use  continuity 
of  df/d€  to  evaluate  q,  df/d?  was  evaluated  for  € < q and 
€ > q , and  the  two  expressions  set  equal.  This  led  to  a 
complex  value  for  q,  which  is  unacceptable.  As  a result, 

1 
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the  following  procedure  was  used  to  choose  q.  Let 
df(€)/de_^  = df/d€  , for  € > q , and  df(€)/d€_  = df/d€  , 

for  ^ < q • Then  let  A(q)  = df(q)/dc_^  - df(q)/de 
Letting  dh/dq  = 0 results  in  a value  of  q that  minimizes 
the  change  in  df/dC  at  q.  The  value  of  q that  resulted  was 

2kTx. 

Validity  of  Assumptions . It  is  fairly  easy  to  examine 
the  validity  of  the  assumptions  involved  in  the  solution  for 
€ < q . As  was  shown  in  Eq  (2-15).  for  ^ < q . 3f/dG  is 
given  by 


M exp[-xATx]  ^ Y /-g  \’ 

" 2kTx 


n!(n  - h) 


(B-1) 


where  the  validity  of  the  solution  is  tested  at  € = q 
The  speed  with  which  this  converges  is  very  sensitive  to  the 
choice  of  q.  As  an  example,  Table  B-1  gives  the  values  of 

An.  where  An  = n!(n  - i)  ■ 'J  ' , 

and  q = 3kT 


n 

A 

n 

q = kTx 

q = 2kTx 

A 

n 

q = 3] 

1 

1.0 

2.0 

3.0 

2 

o 

00 

.3333 

.75 

3 

.0083 

.0667 

.225 

4 

.00074 

.0119 

.06 

Table  B-1.  Convergence  of  af/3€ 
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As  can  be  seen,  the  choice  of  q = kT^  allows  retention 
of  only  the  first  two  terms  in  the  series  with  excellent  re- 
sults . 


€ = X 


It  is  almost  as  simple  to  examine  the  assumption  for 
L . By  comparing  the  exact  expression  for  the  integral 


in  equation  (2-11)  with  approximate  form  actually  used,  it 
can  be  shown  that  the  error  is  proportional  to  A,  where  A 


is  given  by 


/. 


(“x) 


2 -u, 
u e du 


(B-2) 


1 ^ 


By  evaluating  both  A and  the  similar  expression  in  the 
approximate  form,  and  utilizing  the  assumptions  used  through- 
out the  derivation,  it  can  be  shown  that  the  following  condi- 
tion is  sufficient  for  the  approximation  to  be  valid: 


1 » exp 


^ 

4 4 X 


(B-3) 


This  condition  is  much  more  sensitive  to  than  it  is  to 

■^x 

€ X 

X'  kT~  ~ ^ typical  value,  Eq  (C-3)  leads 

to  1 » .16  for  € = 0 , and  1 » .11  for  € = kT^ 

This  approximation  is  not  as  sound  for  q = kT^  as  is  the 
approximation  for  € « X : however,  increasing  q would 
effect  little  improvement  on  the  approximation  for  € > kT^  , 
and  it  would  seriously  degrade  the  approximation  for  6 < kT 


63 


Since  the  physical  and  mathematical  arguments  have  al- 
ready indicated  that  a choice  of  q ® kT^  is  valid;  and 
since  it  has  ,iust  been  shown  that  q = kT^  results  in  valid 
results,  that  value  was  used. 

Three-Level  Model 

Since  the  derivation  for  the  three-level  model  is  exact- 
ly analogous  to  the  two-level  derivation,  the  same  choice  of 
q would  be  , valid,  except  that  there  are  two  different  values 
for  : kT^j^  and  Regardless  of  the  choice  of  q,  f(e) 

should  reduce  to  the  two-level  case  for  and  = 0 . 

Under  these  circumstances,  therefore,  q should  reduce  to  kT^ . 
In  addition,  note  that  f(e)  for  € < q is  affected  strongly 
by  the  behavior  of  f{€)  at  + q and 

Xg  < e < X2  + q (see  Fig  (B-2)).  Since  f(Xj^)  » f(X2)  , 

it  is  reasonable  to  expect  that  f(€)  at  € = 0 is  affected 
most  strongly  by  f(€)  at  € * X^  . For  this  reason  the 
choice  of  q = kT^^^  was  made. 
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A. 


Figure  B-2.  Comparism  of  f(£)  for  € 
and  € 


! 

f 
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1 


NGB  is  a Fortran  computer  program  that  numerically 
solves  the  Boltzmann  equation  for  conditions  typical  of  elec- 
tron-beam sustained  discharge  lasers.  It  can  include  effects 
due  to  field-driven  ionization,  electron-beam-driven  ioniza- 
tion, excitation  and  de-excitation,  heating  by  a constant 
electric  field,  electron-electron  interactions,  elastic  col- 
lisions with  the  background  species,  and  disassociative  at- 
tachment due  to  reactions  of  the  form  e + Cl^  — ^ Cl  + Cl" 

It  integrates  the  Boltzmann  equation  to  find  a time-dependent 
f(€).  It  continues  to  output  f(€)  as  a function  of  time  until 
f(€)  is  constant  from  one  time  step  to  the  next.  Note  that 
since  f(€)  represents  the  fraction  of  the  total  electron  den- 
sity at  a particular  energy,  it  is  possible  to  reach  a quasi- 
steady solution  where  f(€)  is  constant,  but  n^  is  increasing. 

The  integrator  used  in  NGB  is  the  IMSL-designed  DVOGER. 
DVOGER  is  a variable-order  Gear-type  integrator;  in  NGB 
DVOGER  will  go  as  high  as  sixth  order. 

For  a fuller  explanation  of  NGB,  see  the  AFIT  Tech  Re- 
port to  be  published  in  the  fall  of  1978,  by  Seward,  Hastings, 
Ercoline,  and  Gist. 

Particular  Application 

Program  Modification.  NGB  was  modified  slightly  for 
some  portions  of  this  thesis.  The  standard  NGB  calculates 
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the  heating  due  to  an  externally  applied  field  by  calculating 
the  average  collision  time  for  an  electron  of  a particular 
energy  with  all  the  background  species,  and  then  calculating 
the  current  produced  by  those  electrons  under  the  influence 
of  the  external  field.  This  produces  an  energy-dependent  €. 
Since  the  analytic  model  assumed  a constant  €,  NGB  was  mod- 
ified  to  return  a constant  6 = 2 x 10  eV/sec  for  all 
solutions . 

Input  Parameters . NGB  allows  the  user  to  select  which 
effects  he  wishes  to  include.  For  this  thesis  the  input  pro- 
gram switches  were  set  to  include  only  the  effects  of  heating 
(subroutine  HEAT)  and  non-ionizing  inelastic  and  super- 
elastic collisions  (subroutine  XCITE) . An  initial  and  min- 

-14 

imum  timestep  of  1 x 10  seconds  was  used.  The  internal 
error  parameter  for  DVOGER  was  set  at  .001.  The  maximum 
allowed  convergence  difference,  which  defines  the  conver- 
gence condition,  was  set  to  1 x 10 
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Appendix  D 
Programs 

General 

The  following  programs  allow  calculation  of  the  values 
of  the  analytic  two  or  three  level  solutions  on  a Texas  In- 
strument Programmable  59  calculator.  The  program  for  the 
two-level  solution  requires  177  steps;  the  program  for  the 
three-level  requires  477  steps.  Each  program  is  segmented, 
and  could  be  run  a segment  at  a time  on  a TI-58  calculator; 
however,  since  the  TI-58  does  not  allow  automatic  program 
readin,  this  would  be  a time  consuming  procedure. 

Two-Level  Program 

Program  Listing.  Table  D-1  gives  the  program  listing 
for  the  two-level  program. 


step 

156 


160 


170 


176 


Key  Code Operation Comment 


66 

55 

43 

PAUSE 

1 

RCL 

Calculate 

00 

00 

f(€)/N,  for 

95 

= 

e < X 

34 

7x 

94 

V- 

85 

+ 

01 

1 

95 

= 

65 

X 

43 

RCL 

04 

04 

95 

= 

22 

INV 

23 

LNX 

55 

4 

43 

RCL 

05 

05 

95 

= 

91 

r/s  — 

Table  D-1.  Program  for  Two -Level  Model 


Program  Use . First,  X is  stored  in  register  00,  kT^  in 

register  01,  d€  in  register  02,  and  an  initial  value  of  € in 

03.  Pressing  key  A causes  the  calculator  to  evaluate  N, 

storing  it  in  register  05,  and  to  evaluate  the  constant  ex- 

X 


pressions  1 + 


kT 


12 

12  VkT 


in  equation  (2-18)  and 


store  it  in  register  06.  Pressing  key  B causes  the  calcula- 
tor to  increment  € by  d€ , display  the  new  € momentarily,  and 
then  calculate  and  display  f(C),  using  the  expression  valid 
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for  6 < kT 


Fi'essing  key  C similarly  calculates  f(€) 


for  kT^  < e < X : D calculates  f(e)  for  X < € . Note  that 
the  program  does  not  check  to  insure  that  the  value  of  € and 
the  form  of  f(€)  are  compatible. 


Three-Level  Program 

Program  Listing.  Table  D-2  gives  the  program  listing 
for  the  three-level  program. 


Ste 


Key  Code 


Operation 


Comment 


203 


210 


r 

I 220 


230 


54 

) 

42 

STO 

23 

OQ 

02 

65 

X 

Calculate  N; 

53 

( 

store  in  R 09 

43 

RCL 

(here  until 

02 

02 

step  364) 

34 

7x 

65 

X 

43 

RCL 

21 

21 

85 

+ 

02 

2 

65 

X 

43 

RCL 

22 

22 

55 

■ 

• 

43 

RCL 

02 

02 

34 

7x 

85 

+ 

02 

2 

65 

X 

43 

RCL 

23 

23 

55 

43 

RCL 

02 

02 

71 

SBR 

45 

X 

y 

75 

- 

43 

RCL 

02 

02 

65 

X 

80 


step Key  Code 


Operation 


Comment 


378 

55 

4- 

43 

RCL 

380 

11 

11 

54 

) 

71 

SBR 

45 

X 

y 

75 

- 

43 

RCL 

01 

01 

55 

4- 

02 

2 

55 

4- 

390 

43 

RCL 

11 

11 

85 

+ 

01 

1 

54 

) 

85 

+ 

43 

RCL 

40 

40 

55 

4- 

02 

2 

400 

65 

X 

43 

RCL 

02 

02 

71 

SBR 

^5 

X 

y 

65 

X 

53 

( 

43 

RCL 

11 

11 

34 

7x 

410 

55 

4- 

02 

2 

55 

4- 

43 

RCL 

' 

85 

86 


Step 

Key  Code 

Operation  Comment 

449 

08 

08  — 

450 

43 

RCL- 

30 

30 

Calculate 

65 

X 

constant  factor 

43 

RCL 

in  Eq  (2-21), 

01 

01 

for 

71 

SBR 

V 

“xl  ^ « •=  “x2 

45 

A 

y 

store  in  R 07 

55 

•f 

43 

RCL 

11 

11 

460 

85 

+ 

43 

RCL 

40 

40 

65 

X 

43 

RCL 

02 

02 

71 

SBR 

45 

X 

y 

55 

43 

RCL 

470 

21 

21 

54 

) 

42 

STO 

07 

07 

43 

RCL  — 

09 

09 

Display  N 

476 

91 

r/s- 

Table  D- 

Program  Use . 


2.  Program  for  Three-Level  Model 
First,  initial  parameters  are  load 


c , 


Table  D-3. 
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AD-A064  213 


UNCLASSIFIED 


AIR  FORCE  INST  OF  TECH  WRI6HT-PATTERS0N  AFB  OHIO  SCH— ETC  F/6  20/5 
ANALYTIC  electron  ENER6Y  DISTRIBUTION  FUNCTIONS  FOR  ELECTRIC  OI--ETC(U) 
DEC  78  T E 6IST 


AFIT/6EP/PH/78D-4 


NL 


END 


4-79 


Register 


Content 


€ (usually  0 initially) 


^ - x^- 


Rl  + Rg 


Table  D-3*  Initial  Parameters 


Then  key  A'  is  pressed.  The  calculator  evaluates  additional 
parameters,  stores  them  as  shovm  in  Table  D-4,  and  displays 


Register 


Contents 


- x^- 


T R 

^2 

j kT;n^  * * ^2 


Constant  term  in  Equation  (2-23), 


€ < kT 


Term  starts 


and  ends 


F2LLiil/'i2 

L 
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Contents 

(1cT^2)3 


Table  D-4.  Post -Initialization  Parameters 


From  this  point,  pressing  A,  B,  C,  or  D causes  the  calculator 
to  increment  € by  df . and  calculate  f(6),  using  the  expres- 
sion for  € < < € < X2  1 or 

X2  < ^ I respectively,  depending  on  which  key  is  pressed. 
Again,  the  program  does  not  test  to  see  that  6 and  the  expres- 
sion for  f(€)  are  compatible. 


I 


i 

i 
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